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We obtain prompt di-photon cross-section in proton-nucleus collisions in Hamiltonian light-cone 
approach within a hybrid approximation, treating the projectile proton in the parton model and the 
target nucleus in the Color-Glass-Condensate approach. We study in details the di-photon correla¬ 
tions in quark-nucleus and proton-nucleus collisions at RHIC and the LHC. We show that the single 
fragmentation di-photon produces the away side correlations peak, and the double fragmentation 
component of prompt di-photon is responsible for the near-side peak, and the long-range in rapidity 
near-side azimuthal collimation, the so-called ’’ridge” structure. We study the transverse momen¬ 
tum, density and energy dependence of the di-photon ridge and show that it strongly depends on 
the kinematics and saturation dynamics. We show that while di-photon ridge exists at the LHC in 
quark-nucleus collisions, the effect disappears in proton-nucleus collisions at the LHC. At RHIC the 
ridge-type structure persists at low transverse momenta of di-photon even in proton-nucleus colli¬ 
sions. We argue that di-photon correlation measurments in p+A collisions can help to discriminate 
among models and understand the true origin of the observed di-hadron ridge in p-l-A collisions. We 
also show that in addition to the ridge structure, prompt di-photon correlation also exhibits some 
distinct novel features, including the emergence of away side double-peak structure at intermediate 
transverse momenta. 


I. INTRODUCTION 

Recent experimental observations of two hadrons correlations in the relative azimuthal angle A(j) = 4>i — 4)2 and 
in the pseudorapidity separation Arj = 771 — 772 in high-multiplicity proton-proton (p+p) [H and proton(deuteron)- 
nucleus (p+A) [2 t 3| collisions, show a great deal of similarity to that measured in semi peripheral nucleus-nucleus 
collisions In particular, the discovery of the so-called ridge phenomenon, namely the near-side (A^ Ri 0) di-hadron 
correlations which extend to a large pseudorapidiW separation A 77 , in high-multiplicity events selection in both p+p 
and p(d)+A collisions at the LHC and RHIC [ll-[3|; triggered an on-going debate about the underlying dynamics of 
high-multiplicity events. A similar ridge-type structure has also been observed in heavy-ion collisions at RHIC and 
the LHC and is understood as a phenomenon related to hydrodynamical behavior of the quark-gluon-plasma, see e.g 
Refs, d, 13 . The hydrodynamic approach, based on the final-state interaction provides good description [n|, [I2| of 
the high multiplicity p+p data as well. Nevertheless, no compelling argument has yet been given that a small system 
like the one produced in p+p collisions, an order of magnitude smaller than heavy-ion collisions, should exhibit a 
hydro type behavior. On the other hand, the initial-state Color-Glass-Condensate (CGC) approach, the effective field 
theory of low-x partons in the hadronic or nuclear wavefunctions, also provides a qualitatively good description of the 
same data It is thus fair to say that the true physical origin of the ridge in p+p and p+A collisions is still 

unknown, see also Refs. [23 - [3^ . 

The mechanism of prompt di-photon (photons produced in hard scatterings, not from hadron decay) production 
[ 3 ^ in the CGG/saturation framework is quite different than that of di-hadron production [s^. In the CGG picture, 
hadrons are produced in two stages: first soft gluons (or gluon mini-jets) are put on-shell from the projectile wave 
function by directly scattering on a saturated target, and subsequently gluon mini-jets decay to hadrons via a final- 
state hadronization process. The first stage is theoretically under control while the last stage of hadronization can be 
only treated phenomenologically in the small-x framework using the gluon fragmentation functions [ 3 ^ . In contrast to 
gluons, prompt photons do not scatter on the target gluon field, but rather decohere from the projectile wavefunction 
due to scattering of the quarks. Moreover, prompt photons are free from the final-state hadronization processes, and 
possible initial state-final state interference effects. Therefore, prompt photons can be a powerful probe of initial- 
state effects, and in particular the prompt di-photon correlations can provide vital information about the intrinsic 
correlations of partons in the hadronic and nuclear wavefunctions. Such measurments can also shed light on the origin 
of the observed ridge phenomenon in di-hadron production and help us understand whether it is the initial or final 
state effects that play dominant role in formation of the ridge collimation in p+p and p+A collisions. 

Note also that unlike description of di-hadron correlations [sj, which requires the knowledge of correlators of a 
higher number of Wilson lines, the di-photon production cross section (at least in leading order in as) depends only on 
the dipole amplitude [3^ . which is the best understood observable in terms of the high energy evolution. Therefore, 
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within the CGC framework, the theoretical understanding of observables necessary to describe di-photon production 
is more robust compared to di-hadron production. 

The main aim of this paper is to investigate prompt di-photon azimuthal angular correlations in p-|-p and p-l-A 
collisions at RHIC and the LHC. We concentrate on relatively forward rapidities and employ the hybrid approximation 
within the CGC approach. In Ref. 331, we calculated di-photon cross-section in the hybrid approximation at the 


leading-order in the CGC framework employing the soft approximation. Since in the CGC approach the scattered 
quark does not change its longitudinal momentum while propagating through the target, the soft approximation 
amounts to the assumption that the transverse momenta of the produced photons are smaller that the typical trans¬ 
verse momentum of exchanges with the target. Although the soft approximation greatly simplifies calculations, it is 
not applicable in a significant range of interesting kinematics in high-energy collisions. Here, for the first time, we 
provide the full calculation of prompt di-photon production in p-l-A collisions at leading-order in the CGC framework 
without resorting to any further approximation. 

We show that prompt di-photon in high-energy collisions exhibit ridge-type structure in quark-nucleus collisions 
albeit the physics of these correlations is significantly different than that of hadronic correlations discussed previously 
in the CGC framework [Tdl - Eii . In particular the correlations decrease for moderate rapidity separation between 
the two photons 3 units). We investigate the energy, rapidity, transverse momentum, and density dependence 
of prompt di-photon correlations and show that the di-photon correlations have many novel features. We attribute 
the intrinsic azimuthal collimation (the ridge) in prompt di-photon to the double-fragmentation contribution at 
relatively high transverse momentum while the away-side peak to single fragmentation contribution coming mainly 
from relatively low transverse momentum. We also show that the degree of manifestation of these correlations in the 
final state is strongly influenced by the gluon saturation dynamics. 

This paper is organized as follows: In Sec. H we introduce the main formalism for calculating the cross-section 
of inclusive prompt di-photon production in high-energy proton-proton and proton-nucleus collisions in the CGC 
framework and obtain the cross-section within the Hamiltonian light-cone approch at leading-order. In Sec. HI 
we present the results of numerical calculations for correlations of prompt di-photon production in p-l-p and p-(-A 
collisions at RHIC and the LHC. We summarize our main results in Sec. IV. The details of the calculation are given 
in the Appendix. 


II. INCLUSIVE PROMPT DI-PHOTON PRODUCTION IN THE COLOR-GLASS-CONDENSATE 

APPROACH 

In the CGC approach the quantum corrections enhanced by large lo gar ithms of Xjx are systematically re-summed 
incorporating high gluon density effects at low x and for large nuclei |3fil439l | . The inclusive prompt di-photon pro¬ 
duction h + A ^ - 1-72 -I- V in high-energy dilute-dense scatterings was recently calculated in the CGC approach 

in the soft approximation where a dilute projectile hadron (h) interacts coherently with a dense target A and 
produces two prompt photons 71 and 72 . In the leading order approximation, at forward rapidity, a valence quark of 
the projectile hadron emits two photons via Bremsstrahlung and the produced di-photon-|-jet is then put on shell by 
interacting coherently over the whole longitudinal extent of the target, see Fig. [1] 

Here we derive the expression for the prompt di-photon production cross section off a valence quark in the wave 
function formalism without resorting to the soft approximation. The relevant part of the light front Hamiltonian 
which contains the photon-quark interaction is 
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Here {a, s,i} are color, spin and transverse Lorentz indices respectively, {p,q) are transverse momenta and 
and Oi are the quark, antiquark and photon annihilation operators respectively. The terms we omitted here only 
contribute to virtual corrections for the one quark state and are unimportant in the leading order approximation. 
The one quark state dressed by the photons to order has a general form 
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FIG. 1: Typical CGC diagrams contributing to the semi-inclusive prompt di-photon-quark production at leading-order in quark- 
nucleus collisions. The crossed blob denotes the interaction of the projectile quark to all orders with the strong background 
field of the target nucleus. 


Here the momenta have both transverse and longitudinal component, and momenta qi refer to the momenta of photons 
in the wave function. The coefficient functions B and C will be determined from perturbation theory. The function 
B is of order g, while C is of order g^. The outgoing state after scattering on the target is 


\p,s)out = M{1) \p + l,s)+ j^B{p- q,s,s';q,i)\p + l - q,s';q,i) 

i,s' 
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(3) 


where M(l) is the probability amplitude to transfer momentum I from the target to the quark. The momentum 
transfer from the target I is assumed to be purely transverse. As usual in this formalism we rewrite the outgoing state 
in the dressed quark basis 
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(4) 

In this expression the one quark, one gluon state is still bare. To express it in terms of the dressed quark we need to 
use again Eq. @, but only to order g this time, 

\p + l - qi,s';qi,i)D = \p + I - qi, s'-,qi,i) + / ^ B{p + I - qi - q2, s', s”; q2, j)\p + I - qi - q2, s'’]qi,i, q2, j). ( 5 ) 

The dressing corrections to the two photon state do not have to be taken into account, since they lead to corrections 
of at least a power g^ in the outgoing quark wave function. Thus the final expression for the outgoing state in terms 
of the ’’dressed” eigenstate of the Hamiltonian to order g"^ is 
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Let us define for convenience the amplitude 
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The inclusive two photon cross section is (up to normalization) given by the following expectation value in the outgoing 
state 
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Additionally we average over the spin and the color index of the incoming quark. The result is 
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where N{1) = M*{l)M{l) is the dipole scattering amplitude (see Sec. III). Our next step is to calculate the functions 
B and C. We use the general second order perturbative expressions for perturbation V. 
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where we have omitted the terms involving diagonal matrix elements of V, which in our case vanish. This gives the 
following explicit expressions 
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In order to simplify the equations we rewrite these equation in terms of the energy ratios Q defined as 

q+ = (p+; q+ = (3P+; q+ = ^2P+, 

The result is 
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After a simple algebra in the denominators we find 


2( 1 

B{p-q,s,s',q,i) = — ^(Cp, -(?,) 


2 -C 




(15) 


( 16 ) 



5 


C{p-qi - q 2 ,s,s,qi,i,q 2 ,j) = 


' 4 C 1 C 2 Cl (CiPi-gii) 

(p+)2 1 - (1 [C,ip-qiY 


2-Cl 
Cl 




(17) 


(C 2 P - <?2)fc + (C 192 - C2<?l)fc 


C 2 (ClP- 91)^ + Cl(C 2 P- 92)^ - (C291 - Ci 92 )^ 
We can also obtain the following useful identities, 


2 - 2Ci - C 2 
C 2 


Skj^s" s' s' 


2 C 

B{-q,s,s',q,i) - B{1 - q,s,s',q,i) = -gB — 


9i C^i qj 


(CZ - qY 


2-C 




(18) 


[B{-qi,s,s',qi,i) - B{1 - qi,s,s',qi,i)] B{1 - qi - 92 , s', s”; 92 , j) = 


! V 4 C 1 C 2 

P+ 

2-Cl 


9u 

r,2 


Cih — qii 


(C2^ — 92)fe + (Cl 92 — C29l)fe 


9i (Ci^ 9i) J 

2(1 - Cl) - C 2 


Cl 




C 2 




ICkjC 


(19) 


Ci-qi - q2,s,s',qi,i,q2,j) - C{1 - qi - q2,s,s',qi,i,q2,j) = g 
qii 92fe — (Cl 92 — C29i)fe [Cih — qii) 


2 v'4CiC 2 Cl 


P+ 1 - Cl 
{C 2 I — 92)fc + (Cl 92 — C29i)fc 


qi C 29 ? + Ci9i - (C 291 - Ci92)^ [(Ci^ - 91 )^] [C 2 (Ci^ - 91 )^ + Ci(C 2 ^ - 92 )^ - (C 291 - Ci92)^] 


2 -Cl 


Cl 


^li^ss" ss" 


2(1 - Cl) - C 2 
C 2 


^kj^s" s' 


( 20 ) 


Finally we can write 


F(0,Z,9i,92,u,s,s') = g 


2 V4C1C2 


2-Cl 
Cl 


^li^ss" ss" 


2(1-Ci)-C2 

C 2 


^kj^s" s' s" s' 


[ Cl 

qii 92fc — (Ci 92 — C29l)fc 

(Ci^; — qu) {C 2 I — q2)k + (Ci 92 — C29i)fe 

li-Ci 

9 ? [C 29 ? + Ci9i - (C 291 - Ci92)^] 

[(Ci^ - 9i)^] [C 2 (Ci^ - 9i)^ + Ci(C 2 ^ - 92 )^ - (C 291 - Cl 92 )^] 


qu Cik — qu 

{C 2 I - 92)fc + (Cl92 - C 291 ) 

A 

. qf iCii-qiY. 


{C 2 I - 92 ) + (Ci92 - C 291 ) 

A 


( 21 ) 


Our goal now is to calculate the production cross section using Eq. This involves symmetrizing the amplitude F 
in Eq. (1211) and squaring it. The algebra is fairly lengthy, and we present the derivation in the Appendix. The final 
result is 
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where the functions F^p and F^p are given in Eqs. (I36I37I38I39|) in the Appendix. Note that in the final expersions 
for Fap and F^p in Eqs. (I36I37I38I39|) all momenta q,ki,k 2 ,l are 2-dimensional vectors in the transverse plane. For 
clarity, in Eq. dm) and in the following, transverse momentum vector is denoted by a bold notation with a subscript 
T. 

In the next section we discuss our numerical analysis of Eq. ([2^ . The simple case of the eikonal (soft) approximation 
where the analytical expersions can be significantly simplified, is discussed in the Appendix. 
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III. NUMERICAL RESULTS AND DISCUSSION 


The dipole amplitude N{It) that enters Eq. (|2^ defends on energy through the longitudinal momentum variable 


Xgilr) ! which is defined as follows (see Appendix in [33|), 
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where the parameter Xq is the ratio of the incoming quark to the projectile nucleon energy with ^/s being the nucleon- 
nucleon center-of-mass energy. The parameter Xg is the fraction of longitudinal momentum of the target carried by the 
gluon which is absorbed by the scattered quark. As discussed in the last paper in [s^ this value of Xg is appropriate if 
one assumes that all the momentum from the target is transferred to the quark in a single scattering event. Otherwise 
one should use a different value of Xg consistent with multiple scattering off a dense target. However this issue is only 
important at NLO accuracy in as, and in the present leading order calculation we follow the time honored practice 
and use the value given in Eq. (1^51) . The dependence of the dipole scattering amplitude on Bjorken Xg is determined 
by the JIMWLK renormalization group equations. 

The function N {lT,Xg) in Eq. (|22)l is related to the correlator N{r,Xg) of two light-like fundamental Wilson lines 
in the background of the color fields of the target nucleus (or proton) through the Fourier transform. In the large Nc 
limit, the coupled JIMWLK equations [s^ are simplified to the Balitsky-Kovchegov (BK) equation [s^, a closed-form 
equation for the rapidity evolution of the dipole amplitude which is presently known to next-to-leading accuracy 
( 40 I Idljl . The running-coupling improved BK equation (rcBK) has the same generic formal form as the leading-log 
BK evolution equation: 
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where r = |r 2 — ri| is the transverse size of dipole and the is evolution kernel with Balitsky‘s prescription (43| 
for the running coupling. Note that inclusive di-photon production cross-section (both direct and fragmentation part) 
depends on the dipole-target amplitude and therefore in principle probes the small-x dynamics. In the following we 
study the quantity, which is related to the above partonic cross section via 
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where fq denotes the parton distributions function (PDF) in the projectile hadron (nucleon). The lower limit of 
integral x™*" in the above is defined by, 
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Denoting the quantity in Eq. (EH) by the two photon differential cross section is a certain abuse of notation, since 
the actual cross section contains also contribution of photons emitted independently from two quarks (antiquarks). 
This extra contribution, although of the same order in aem and as is of little interest to us, since we do not expect 
any correlated emission from independent quarks. It will likely amount to an uncorrelated “pedestal” proportional 
to the product of single photon inclusive cross sections, which will have to be subtracted from the data (if and when 
such data is available) in order to compare with our results. Formally this subtraction can be written as 


(l(jPA^l{ki)l(k2)X ^^pA^'r{ki)'f(k2)X ^^pA^'y(ki)X ^^pA^j(k2)X 
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(27) 


The effective area cr(|fci|, |fc 2 |, ? 7 i, 772 ) is a slowly varying function of momenta (when the momenta are large enough) 
and is of the order of the area of the nucleus (or target proton). However its calculation from first principles requires 
an additional nonperturbative input, and rather than attempting it we suggest to treat it as an adjustable parameter 
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FIG. 2: Left: the Xq (ratio of energies of the projectile quark to nucleon) dependence of prompt di-photon correlations C 2 
obtained by solving the rcBK evolution equation with an initial saturation scale for the target QosA = 5Qos with Qqs = 
0.168 GeV^ corresponding to minimum-bias quark-nucleus collisions. Right: The initial-saturation-scale dependence of prompt 
di-photon correlations C 2 in q-|-p(A) collisions as a function of angle between two photons Acp. The results are obtained at a 
fixed Xq = 0.01 by solving the rcBK evolution equation with different initial saturation scales for the target QqsA = NQQg with 
A’ = 1, 5,15 and Qqs = 0.168 GeV^. We also compare with the non-saturation (No Sat) model, see the text. All curves in both 
panels are obtained at a fixed energy ^/s = 5.02 TeV, at fixed rapidities of photons 771 = 1 r ;2 = 2 with transverse momenta of 
photons taken to be kir ~ k 2 T = 2 GeV. 


in comparing with future data. In the rest of this paper we retain the notations of Eq. (HSl), while keeping the above 
caveat in mind. 

The only external input for the rcBK non-linear evolution Eq. (1241) is the initial condition for the evolution which 
is taken to have the following form motivated by McLerran-Venugopalan (MV) model 


M{r, V = 0 ) = 1 — exp 




(28) 


where the onset of small-x evolution is assumed to be at ccq =0.01, and the infrared scale is taken A = 0.241 GeV (43 |. 
The only free parameters in the above are 7 and the initial saturation scale Qos- The initial-condition of the rcBK 
equation has a non-perturbative nature and its parameters are fixed via a a global fit to proton structure functions 
in DIS in the small-x region [4^ and single inclusive hadron data in minimum-bias p+p collisions at RHIC and the 
LHC [dst - lrij . Note that the current HERA data alone are not enough to uniquely fix the initial condition of the 
rcBK. Different parameter sets of the rcBK equation correspond to different initial saturation scale Qos (as probed 
by quarks). For proton-proton collisions, we take the initial saturation scale of proton ~ 0.168 GeV^ (with the 
corresponding 7 ~ 1.119) which was extracted from a global fit to proton structure functions in DIS in the small-x 
region and single inclusive hadron data in p-l-p collisions at RHIC and the LHC [isl - li^ . For the nucleus case, 
the initial saturation scale of a nucleus Qg^ ~ 5 A 7Qos should be considered as an impact-parameter averaged value 
and it is extracted from the minimum-bias data in deuteron-gold at RHIC and proton-lead collisions at the LHC . 
In the absence of proper impact-parameter dependence of the dipole amplitude, rare events with higher multiplicity 
than the minimum-bias, may be described by using a hi ghe r initial saturation scale incorporating somehow the effects 
of fluctuations and geometry [H, IH, see also Refs. [49l - l^ . 

Let us define the azimuthal correlation of the produced di-photon as [H, [sj , 
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where A(j) is the azimuthal angle between the two produced photons in the plane transverse to the collision axis. 
The azimuthal correlation C 2 is proportional to the relative probability of inclusive di-photon pair production in a 















FIG. 3: Prompt di-photon correlations C 2 in quark-nucleus collisions as a function of angle between two photons A(^ in q-|-Pb 
collisions at the LHC {y/s = 5.02 TeV) with ? 7 i = 1,772 = 2 and Xq = 0.01. The curves are results of the rcBK solution with 
two initial saturation scale, similar to Fig. [2] right panel. A schematic contribution of different typical CGC diagrams to the 
double inclusive prompt photon production at different angle is also shown. 


given kinematics and angle between the photons in the pair, normalized to the total cross-section integrated over 
angle Following the experimental procedure for the zero-yield-at-minimum (ZYAM), we remove the constant 
background Czyam by shifting the minimum at the zero axis. We expect that some of the theoretical uncertainties, 
such as sensitivity to possible AT-factors which effectively incorporates the missing higher order corrections, drop out 
in the correlation defined in Eq. One can equally take the normalization in Eq. as the differential cross- 

section at a fixed reference angle A^c (away from near-side or away side). We checked that this does not alter our 
over-all picture here. One can equally define the correlation via coincidence probability however since the 

correlation defined via Eq. (I29p is free from extra integrals over transverse momenta, it should exhibit the underlying 
dynamics of the correlation in a cleaner way. Note that defining the correlation for prompt di-photon via coincidence 
probability depends on the definition of the trigger single inclusive prompt photon whether it is taken to be direct or 
fragmentation prompt photon or include both contributions. The correlation definition given in Eq. (ESI) is free from 
this ambiguity. 

For consistency of our leading order approach, we should use the leading order PDF’s. However the ratio Eq. ESD 
depends on PDE’s only very weekly. We checked that the leading-order (LO) and next-to-leading-order (NLO) PDFs 
give similar results for the di-photon cross-section, with less than 10% discrepancy (which can be absorbed into the 
AT-factor). In the correlation C 2 defined via Eq. (E^ the effect of varying PDE’s between LO and NLO is practically 
negligible within less than 2% uncertainty. The results we present in the following were obtained using the NLO 
MSTW PDEs [ 53 . Eollowing the conventional pQCD approach, we take the hard scale /r/ = {kiT + kiT)l2 in the 
parton distribution in Eq. (1251) . For numerical computation, we focus here on low transverse momenta of the produced 
photon pair at the LHC at forward rapidities. Note that this kinematics is mostly relevant for probing saturation 
effects and as we will demonstrate here, it is in fact the kinematics in which the ridge structure is by far the most 
pronounced. 

We first focus on the prompt di-photon correlations in quark-nucleus (q-l-A) collisions. In Fig. EJ left panel, we 
show the di-photon correlations C-i defined via Eq. ESI) in quark-nucleus collisions at the LHC energy 5.02 TeV at 
fixed rapidities 771 = 1 and 772 = 2 for different values of Xq. We recall that the parameter Xq is the ratio of the 
incoming quark to the projectile proton energy. All curves in Eig. [2] left panel were obtained by solving the rcBK 
evolution equation with an initial saturation scale for the target QqsT = ^Qos Qg^, = 0.168 GeV^ corresponding to 
minimum-bias collisions. It is seen that the near-side and away-side correlations (the peaks) is enhanced by decreasing 
Xq. For a large Xq, all correlations both at near-side and away-side disappear, albeit this effect is more pronounced for 
the near-side correlations. Note that at fixed rapidity, energy and transverse momentum, a large Xq corresponds to 
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FIG. 4: The prompt di-photon correlations in q-l-A collisions as a function of angle between two photons A(j) obtained by 
imposing different lower cut in the total transfer momentum to target (or total transverse momentum of the produced quark- 
diphoton) denoted by The results are obtained within the rcBK saturation model with the initial saturation scale 

QosA ~ ISQos and non-saturation model with leading-order perturbative unintegrated gluon density. All curves are obtained 
at ^/s = 5.02 TeV with Xq — 0.01, and 771 = 1, ?72 = 2 and transverse momenta of photons fcir = kir = 2 GeV. 


small Cl and C 2 , see Eq. (1231) . Our full numerical calculation for C 2 shown in Fig. [2] is consistent with our theoretical 
conclusion about the absence of correlations at small Q discussed in the Appendix. 

In Fig. [21 right panel, we show the initial saturation scale dependence of C 2 at a fixed Xq = 0.01 in quark-proton and 
quark-nucleus collisions at 5.02 TeV with rapidities of the produced photon rji = 1 and 772 = 2. Different curves were 
obtained by solving the rcBK evolution equation with different initial saturation scales for the target = NQos 

with N = 1,5,15, and Qq^ = 0.168 GeV^ being the initial-saturation scale of the proton. We also compare with 
the non-saturation model (labeled by ”No Sat” in the plot), by taking the unintegrated gluon distribution at the 
leading-order pQCD, N{It, Xg) oc l/^|n in the prompt di-photon cross-section. It is seen that increasing the saturation 
scale in the system increases the near-side correlations leading to near-side collimation, while it decreases away- 
side correlation leading to the away-side decorrelation. Note that if momentum transfer from the target is small 
enough, the photon collinear to the outgoing quark, and the photon emerging from the initial photon-quark state 
have opposite transverse momenta, leading to back-to-back correlation at = tt due to the single fragmentation 
di-photon contribution. Now increasing the exchanges to the target by increasing the saturation scale, unbalances 
the back-to-back correlations. The double fragmentation and direct di-photon contributions (see Ref. [s^) on the 
other hand, are restricted to kinematics where the transverse momentum of the outgoing quark jet is relatively large. 
One therefore does not expect significant back-to-back correlation in the double fragmentation and direct di-photon 
contributions. By increasing the saturation scale, the main contributions of integrand is shifted to higher transverse 
momentum, leading to the enhancement of double-fragmentation contributions at the near-side Ac/) = 0. The typical 
CGC diagrams which contribute at different A^ are also shown in Fig. [31 The curves in Fig. |3| are the same as in 
Fig. [2] right panel. 

In Fig. m we show the effect of lower cut on the total momentum transfer to target It in Eq. (1331) (in the cross- 
section at partonic level) on the correlations C2 within the saturation and non-saturation models. For the saturation 
model, we use the rcBK solution with initial saturation scale Qg^^ = 15Q§g corresponding to a rare high-multiplicity 
event. All curves in Fig. 0|are obtained at a fixed energy ^/s = 5.02 TeV, at fixed rapidities of photons 771 = 1 and 
772 = 2 with transverse momenta of photons taken to be kiT — k2T = 2 GeV. First it is seen that both saturation 
and non-saturation models produce similar correlations with a lower-cut It > = 4 GeV. This is clear since the 

saturation scale of the system in the chosen kinematics region is smaller than and therefore saturation effects 

should not be important for It > It^'^- Consequently both the saturation and the non-saturation models give similar 
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FIG. 5: Left: the rapidity dependence of the ratio of prompt di-photon differential cross-section in proton-nucleus collisions 
as a function of angle between two photons Acj), normalized to integrated cross-section over the angle. The curves are obtained 
for p-|-Pb collisions by the rcBK solutions with the initial saturation scale QqsA = 15Qos Right: the initial-saturation-scale 
dependence of prompt di-photon correlations in p-|-p(A) collisions. The results are obtained by using the rcBK evolution 
equation solutions with different initial saturation scales for the target QqsA = a-nd = 1,5,15. We also compare with 

the non-saturation (No Sat) model. All curves are obtained at y/s = 5.02 TeV, and r]i = 1, t ]2 = Arj + t]i and transverse 
momenta of photons fcir = fciT = 2 GeV. 


results. It is clearly seen in Fig. 0] that increasing (from 0.2 GeV to 4 GeV) enhances the near-side correlations 
while decreases the away-side correlations. As we expect, the effect of changes on total momentum transfer to target 
(or changes on the lower-cut 1 ™*") on the correlations is qualitatively similar to the effect of variation of the saturation 
scale Qs of the system. 

Recall that as shown above, the di-photon near-side collimation appears in quark-nucleus collisions only at small 
value of Xq. In order to compute the correlations in proton-nucleus collisions, we should convolute the partonic cross- 
section with the PDFs of quarks and anti-quarks, and perform the integral over Xq, see Eq. (ESI) . Note that although 
the partonic cross-section explicitly depends on l/x^, in proton-nucleus collisions this will be modified due to Xq 
dependence of the PDFs. Moreover, the convolution integral in Eq. (1251) has a lower limit on Xq via Eq. (l26l) which 
depends on the energy/rapidity and transverse momentum of the produced photon. Note also that the parameter Xq 
also appears in the definition of Xg via Eq. (l23ll and hence the PDF convolution indirectly affects the rcBK evolution 
solutions of target nucleus in a non-trivial way (mainly due to energy-momentum conservation constraint). Therefore, 
a priori, it is not obvious how the integration over Xq will affect the di-photon correlations in p-l-A collisions. 

In Fig. [5] left panel, we show the rapidity dependence of the ratio of prompt di-photon differential cross-section in 
proton-nucleus collisions as a function of angle between two photons Ac/), normalized to integrated cross-section over 
the angle. The curves are obtained for p-l-A collisions at the LHG energy 5.02 TeV with fixed transverse momenta 
of photons fciT = fciT = 2 GeV. The initial saturation scale of the target was taken = ^Qos ^ = 15, 
corresponding to an event selection with high charged hadron multiplicity. Recall that within the same multiplicity 
event selection (with N = 15), the di-hadron ridge type structure is well pronounced [l^. On the other hand we see 
here that within the same multiplicity event selection, the prompt di-photon near-side collimation almost disappears 
in p+A collisions at the LHC. In order to examine if the di-photon near-side correlations survive at any multiplicity 
event selections, in Fig.[5]right panel, we show the initial-saturation-scale dependence of prompt di-photon correlations 
defined via C 2 in p-l-p(A) collisions at 5.02 TeV and 771 = 1, 772 = 3 with fixed transverse momenta of photons 
kiT = kiT = 2 GeV. The effect of size and density of the target is simulated by taking different initial saturation scale 
for the rcBK evolution equation. Namely we use the different solutions of the rcBK evolution equation with different 
initial saturation scales for the target Qqsa = ^Qos ^ — 1; 5,15. We recall that an initial saturation scale Qos 
describes a proton target in minimum-bias electron-proton or proton-proton collisions , while an initial saturation 
scale of about 5 -L 7Qos describes a heavy nuclear target in minimum-bias proton-nucleus collisions (dsl - l^ . In order 
to examine the effect of saturation more clearly, we also show the results obtained by taking the dipole amplitude 
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FIG. 6: The rapidity (left) and transverse momentum (right) dependence of the prompt di-photon correlations in proton- 
nucleus collisions at 5.02 TeV as a function of angle between two photons Acj) obtained within the non-saturation model with 
a lower cutoff on the total transverse momentum of the produced quark-diphoton taken to be Zy = 4 GeV. 


at the leading-order pQCD (labeled by ”No Sat” in the plot). It is seen from Fig. [Sj while increasing the saturation 
scale suppresses the away-side correlations, it does not affect the near-side correlation at the LHC energy 5.02 TeV 
in p-|-A collisions. 

In order to examine the effect of a lower cut on the total transverse momentum Z™*", similar to calculations shown 
in Fig. m we imposed a lower cut Z™™ = 4 GeV, and we computed the correlation C 2 in the presence of such a cut in 
p-|-A collisions at the LHC energy 5.02 TeV. In Fig. [6l we show the rapidity dependence (left panel) and transverse 
momenta of di-photon dependence (right panel) of the prompt di-photon correlation C 2 in proton-nucleus collisions 
within the non-saturation model with a lower cutoff. It is again seen that similar to the case of quark-nucleus collisions, 
imposing a lower cut Z™“, enhances the near-side correlations. Note that in principle double photon fragmentation 
should not contribute when the rapidity of the produced photons are different^. However, the remnant of the effect of 
double fragmentation may still persist up to a couple of units in the rapidity difference (depending on the kinematics 
and energy). This is shown in Fig. |6]left panel where we show C 2 as a function of Ayy. It is seen that the correlations 
both in near-side and away-side in p-l-A collisions at the LHC disappears for approximately At; > 3. Therefore, 
these correlations are relatively short-range in nature and may not be considered as a genuine ridge compared to the 
di-hadron ridge. In principle, in order to enhance these correlations, one can choose a kinematics region more sensitive 
to the double fragmentation contribution. This can be achieved by triggering quarks (or hadrons) with a relative 
large transverse momentum simulating the effect of a large Z™*". The exact value of Z^ depends on kinematics (and 
the acceptance of detector). The systematics of our results and the kinematics that the long-range correlations can 
be probed, will be outlined at the end of this section. 

Note that the rapidity, energy and density dependence of di-photon cross-sections and its evolution via the rcBK 
equation is very non-trivial due to complicated form of Xg in Eq. (1231) that enters in the dipole amplitude in the 
cross-section. One also should bear in mind that generally, the small-x evolution in the rapidity at a fixed energy 
and transverse momentum is not the same as the evolution in energy but at a fixed rapidity. This is partly due to 
the fact that the available phase space for particle production at different energy is different. Having this in mind, 
we next focus on prompt di-photon correlations at RHIC. In Fig. fright panel, we show the transverse momentum 
dependence of the prompt di-photon correlations C 2 in minimum-bias proton-nucleus collisions obtained with the 
rcBK solutions with the initial saturation scale (simulating minimum-bias p-l-A collisions with a heavy 


^ In principle, similar to Ref. [ 3 ^ . one can explicitly extract from the prompt di-photon cross-section, single and double fragmentation 
functions and contributions, see Appendix B. However, this is not needed here as we stay away from the collinear kinematics. 
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FIG. 7: Left: transverse momentum dependence of the prompt di-photon correlations C 2 in minimum-bias proton-nucleus 
collisions as a function of angle between two photons A0 obtained with the rcBK solutions with the initial saturation scale 
Q1)sA ~ 5Qos at the RHIC energy 0.2 TeV with 771 = 1 and 772 = 3. Right: comparison of C 2 at RHIC and the LHC for fixed 
transverse momenta of di-photon at kiT = k\T = 1 GeV with 771 = 1 and 772 = 3.5 in minimum-bias proton-nucleus collisions. 


nuclear target) at the RHIC energy 0.2 TeV with 771 = 1 and 772 = 3. In Fig. [7] left panel, we compare the correlation 
C 2 at RHIC and the LHC at fixed transverse momenta of di-photon kiT = kix = I GeV with rapidities 771 = 1 and 
772 = 3.5 in minimum-bias proton-nucleus collisions. Note that all results in Fig. [7] are obtained with imposing a 
lower cut Z™” ^ Aqcd and the upper-cut ~ 20 GeV in the integral over It at both RHIG and the LHC. It is 
interesting that at RHIC in stark contrast to the LHC, a ridge-like structure appears at low transverse momentum 
of the produced photons 0.2 < ferlGeV] < 3. The near-side collimation at RHIC survives up to rapidity interval 
A 77 Ri 3. Another interesting features of the correlations at RHIC is the emergence of a double-peak structure for 
away-side correlations at low transverse momentum kx of photons. At higher kx it is replaced by a single peak 
structure due to the fragmentation contribution. Note that a double-peak structure for the two-particle structure was 
also reported for the photon-hadron Is^l and dilepton-hadron correlations in a small windows of kinematics in 
the CGC approach. Therefore, it seems the existence of a double-peak structure at away-side to be a universal feature 
of electromagnetic probes. In the present case we understand its emergence in the following way. We recall that in 
order for higher Fock components of the projectile hadron wavefunction to be resolved and photons to be radiated, the 
projectile quark should interact with small-x target by exchanging transverse momentum. This means that when total 
momentum transfer lx is zero, the production rate of di-photon should go to zero and off-shell di-photon remains as a 
part of projectile hadron wavefunction. Note however because of convolution with parton distribution functions, the 
local minimum at low lx will not be zero but gets smeared out. On the other hand, the production has a maximum 
at low lx due to single fragmentation contribution. Therefore, a double peak structure appears for the away side 
di-photon correlations as a result of combination of a local minimum and a local maximum for the correlations. At 
higher transverse momentum, the back-to-back fragmentation contribution overwhelms the cross-section and prompt 
di-photon correlations has a single-peak structure. Note that in the LHC kinematics, the typical momentum transfer 
to the target is larger than at RHIC and the double-peak structure becomes single-peak structure for a wide range 
of kinematics. This effect was also seen for the photon-hadron [H, and dilepton-hadron correlations at RHIC 
and the LHC. 

The systematics of the results can perhaps be understood in the following way. First, recall that in our approach 
the photons do not scatter off the hadronic target directly. The final state di-photons arise from two contributions. 
One is the two photon component of the incoming quark state, which is put on shell by virtue of the quark scattering 
out of the initial state, while the other is the photon(s) emitted by the ’’bare” outgoing quark while it ’’reconstitutes” 
the dressed quark wave function, see Fig.[T] The first contribution is a part of the direct di-photon, while the second 
is mostly contained in the fragmentation (single and double fragmentation) part, see Fig. [31 There is no physics in 
the scattering event itself that can generate the same side correlations. Therefore such correlations can only appear 
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due to correlations that exist in the proton wave function or are produced during the (double) fragmentation stage. 

Just like in hadronic case Bose enhancement of photons exists in the incoming quark (proton) wave function [59j. 
However now this enhancement does not extend over large interval of rapidities but is peaked at A 77 ^ 2. Once the 
quark scatters out of the initial wave function putting the two photons on shell, the initial state enhancement due to 
double fragmentation translates directly into the same side correlations in the final state, since the momentum of the 
two photons is not changed by the scattering. In this sense the initial state correlations are even more directly reflected 
in the hnal state di-photon momentum distribution than in a di-hadron, which does absorb transverse momentum 
from the target. 

This above simple reasoning assumes that one can cleanly separate the final state photons which originate from 
the two photon initial state. This is difficult to do in the non-saturation model, the ”No Sat” curves on Figs. (IMl 
in]) or when the total transfer momentum to target is very small ^ ^qcd (in the non-saturation model) in 
Figs. (ITO . The scattering probability N{It) ^ 1/?^ is strongly peaked at low momenta. On the other hand, if the 
momentum transfer to the quark is much smaller that the transverse momenta of the observed photons kx -C kix, k 2 T, 
while propagating through the target the quark changes direction adiabatically, so that no photons with momenta 
kiT ^ k 2 T are produced. The production is therefore dominated by the momentum transfer kx kix, k 2 x- In 
this regime however, there is no clear picture of a well defined contribution of the initial two photon state for two 
reasons. First, the factorization scale for fragmentation functions is itself of the order of the momentum transfer 
fi ^ kx and with all momenta being of the same order some of the would be ’’initial state di-photons” are identified 
as fragmentation photons. Secondly, in this regime there is strong interference between different diagrams in the 
amplitude, and the separation on the probability level is not possible. As the saturation momentum is increased 
Qs ^ kix,k 2 x or > Aqcd, the maximum of the scattering probability moves to higher transverse momentum 
transfer kx ^ Qs- At this point the simple picture of initial state correlations becomes valid, and one indeed observes 
the imprint of the initial correlations in the emitted di-photon distribution. Eventually, when Qs and therefore typical 
momentum transfer kx becomes too large, the quark is kicked out of the proton wave function too violently, so that 
an ever increasing contribution to the direct di-photon comes from the immediate vicinity of the interaction point 
rather than from the incoming wave function. Those photons are not correlated, and so the normalized correlation 
function C 2 decreases. 


IV. CONCLUSION 

In this paper, we investigated prompt di-photon correlations, and in particular the existence of ridge correlations 
for di-photon production within the CGC approach in high-energy q-l-p(A) and p-l-p(A) collisions at RHIC and the 
LHC. Our formulation is only valid for the case of dilute-dense scatterings at forward rapidity. Our main findings are 
as follows: 

A ridge-type structure exists for prompt di-photon production in high-energy quark-nucleus collisions (even in 
minimum-bias collisions) at both the LHC and RHIC. Prompt di-photon ridge at near side A(j) Ri 0 mainly comes 
from the double fragmentation di-photon contribution and the effects of Bose enhancement in the projectile wave 
function, while the away side ridge A(j) Ri tt originates from the single fragmentation contribution which is dominated 
by the usual back-to-back correlations. The ridge structure at near-side in quark-nucleus collisions disappears for 
rapidity differences A 77 ^ 3 . Additionally it exists only in a limited kinematics windows with small value of the 
light-cone variable Xq and within the saturation regime, see Figs. (IMl) . 

For the case of proton-nucleus collision, the near-side correlations are washed away at the LHC 5.02 TeV after 
integration over Xq variable (or convolution with the PDFs) while the ridge-like structure survives in p-fA collisions 
at RHIC 0.2 TeV at rather low transverse momenta of the produced photons even up to Arj Ri 2 ^ 3, see Figs. (I5I7I) . 
We showed that in principle the ridge-like structure at the LHC can re-appear by imposing a lower cut on the total 
transfer momentum to the target by isolation cut techniques, see Fig. [51 

The prompt di-photon near side At)) k, 0 correlations (the ridge) strongly depends on the saturation dynamics and 
the value of typical total transfer momentum to the target (or It is absent when Qs kx, is generated by 

increasing the saturation scale to kx < Qs or kx < small Xq). This is mainly due to the fact that a larger 

saturation scale shifts the main contribution of the integrand of the cross-section to higher transverse momentum 
where the double-fragmentation becomes important leading to enhancement of near-side correlations while it reduces 
the single fragmentation contribution, consequently unbalancing the back-to-back correlations leading to suppressions 
of the way-side correlations, see e.g. Fig. [3| 

The prompt di-photon correlations at forward rapidities in the away side region exhibit single or double peak 
structure depending on the kinematics, namely transverse momenta of the produced photo ns, rapi dity interval between 
two photons and center-of-mass energy. Similar effect was also reported for photon-hadron [53Ll54| and dilepton-hadron 
correlations in high-energy p-fPb collisions. 
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Although di-hadron ridge was observed in p-l-p and p(d)-|-A collisions [IH3 both at the LHC and RHIC, the 
di-photon correlations have not yet been measured at the LHC and RHIC. Our study here shows that the energy, 
density and transverse momenta dependence of di-photon correlations both at near side and away side provide essential 
complementary information to understand the true underlying dynamics of emergence of ridge phenomenon in high- 
energy collisions. It was argued that the di-hadron ridge effect can be understood entirely within the framework of the 
initial-state CGC physics. There are also competing mechanisms based on the final state effects, like Hydrodynamics 
approach [Hill, which provide excellent description of the same data. Our study here pertains to the minimal 
bias p-Pb data, where one certainly does not expect collective phenomena in the final state to play a crucial role. 
Therefore observation or non observation of di photon correlations at RHIC and the LHC (see e.g. Fig. [3) would be 
an important test for the two approaches. 

It is interesting to note that Ref. [l^ required both the projectile proton and the target to be in the saturation 
regime in order to describe the ridge effect in p-fPb collisions. The existence of two saturation scales in the problem 
(for the projectile proton and the target) both in p-l-p and p-fPb collisions, is essential for [I^ in order to find a 
good fit to the observed number of charged hadron tracks in which the ridge appears. The origin of such rare high 
multiplicity events in p-l-p and p-fPb collisions is not yet well understood. In our hybrid approach on the other 
hand the projectile proton is assumed to be in dilute regime and is consequently treated in standard parton model 
approach, while the target is treated as a dense object. Therefore, we have here only one saturation scale and our 
calculation should be applicable to minimal bias events in p-fPb collisions and a more significant fraction of p-fp 
events than those exhibiting di-hadron ridge correlations. In this sense the di-photon correlations is a much cleaner 
probe of initial state effects compared to di-hadrons, as it does not require modelling of rare fluctuations leading to 
high multiplicity events. 

One of the main results of this paper is that the prompt di-photon near-side collimation and the ridge-type structure 
can also be probed in quark-nucleus collisions at the LHC at small value of ccq, albeit within a limited rapidity interval. 
We recall that strictly speaking our formulation is only valid for dilute-dense scatterings (such as proton-nucleus 
and electron-nucleus collisions) at small-x. Therefore, our general outcome here in particular the di-photon ridge 
phenomenon in quark-nucleus collisions can be also tested in a future electron-ion collider [^, . A detailed study 

of di-photon correlations in high-energy electron-nucleus collisions is postponed in future publication. 
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APPENDIX 

In this Appendix we present the details of the calculation of the di-photon production cross-section following the 
formalism introduced in Sec. 11. Symmetrizing Eq. (EH) we obtain: 
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We will calculate the square of the first three lines in Eq. (l30ll and the product of the first three lines by the last three 
lines. The complete result will then be given by adding the symmetrizing term. In the first part of the calculation we 
encounter the expressions of the type 


Ai = —tr 


(2-Ci)(2(I-Ci)-C2), 


X 2 = -tr 


C 1 C 2 

(2-Ci)(2(1-Ci)-C2) 

C 1 C 2 

(2-Ci)(2(I-Ci)-C2) 




C 1 C 2 

(2-C2)(2(1-C2)-Ci) 

C 2 C 1 


Cidj ^ik^k^jldl 
O^ibj f^ik^k^jl^l 


Cidj CikCkCjldl 


6 + i 
S — i 
d + i 
S-i 


2 - Cl ^ , 2(1 - Cl) - C 2 ^ 

CliCjkOk H-T- CilCllOj 


Cl 
2-Cl 


C 2 

, 2(1 - Cl) - C 2 , „ , 

V CiCjkClk \ y. 

. SI (,2 

2 - Cl. . . , 2(1 - Cl) - C 2 _ , 


I Cl 

2(1 - C 2 ) - Cl 


Cl 


C 2 

CiCjkdk ~ CilCldj 

C 2 


where a, b, c, d are various transverse momenta. After some tedious algebra this can be explicitly written as 

Ai = E(12;12)(o-c)(5-(i)-D(12;12)(axc)(5xd), 

A 2 = E(12;21)(a-c)(5-d)-D(12;21)(axc)(5xd). 

Note that the second part of the above expressions can be further simplified by using the following identity, 

(a X c){b X d) = (a ■ b)(c ■ d) — (a ■ d)(b ■ c). 

Let us introduce the following notation 
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[(2-Ci) 2+C?] [(2 - 2Ci - C 2 )" + Cl] 

CICI 

4(2-Ci)(2-2Ci-C2) 

C 1 C 2 

[(2 - Ci)(2(1 - C 2 ) - Cl) + Cl] [(2 - C 2 )( 2(1 - Cl) - C 2 ) + Cl] 

CICI 

4(2 - Cl - C2)^ 

C 1 C 2 


(31) 


(32) 


(33) 


(34) 


(35) 
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We multiply out the different terms 


J^n=/^S(12;12) 

Cl f _ [g2 ~ (Clg2 ~ C2gi)]^ _I_ [(C2^ ~ g2) + (Clg2 ~ C2gl)]^ _ 

(1 - Ci)^ \ql [C2g? + Ci9i - (C 291 - Ci 92)^]^ [(Ci^ - 9i)^] [C2(Ci^ - 9i)^ + Ci(C2^ - 92)2 - (C2gi - Ciq2)'^]'^ 

_ 2[gl • (Cl^ ~ gl)] [(g2 - (Cig2 - C2gi)) • ((C2^ - g2) + (Clg2 - C2gl))] _\ 

9i [C2gi + Cl92 - (C291 - Cl92)^] [(Ci^ - 91)^] [C2(Ci^ - 91)^ + Ci(C2^ - 92)^ - (C291 - Cl92)^] j 

^ _ CE __ 

9?(Ci^ - 9 i)^((C 2/ - 92) + (Cl92 - C291)) 

2(1^ / Ci^ • (Ci^ - 9i) ((Cl92 - C291) - 92) • {{(,2l - 92) + (Cl92 - C291)) 

^ Wi[Ci^ - 9i]^ [C29? + C19I - (C291 - Ci 92)2] {{C2I - 92) + (Ci 92 - C291)) 


_ Ci(^ • gi) _ 

9?[Cl^ - 91 ]^ [C2(Cl^ - 91 )^ + Cl(C2^ - 92 )^ - (C 29 I - Cl92)^] 


A"i = -5"^S(12;12) 

q2 C 1 C 2 (91 X 0 (1 ~ Ci)( 92 X /) + C2(9l X 1) 

" (1 - Cl)" ql [(Ci^ - 91 )"] [C 29 ? + Ci9i - (C 291 - Cl92)2] 

C 1 C 2 (91 X 0 (1 — Ci)(92 X 1) + C2(9l X 1) 

9? (Ci^ - 9i)^ [C 291 + Ci9i - (C 291 - Ci92)^] ((C 2 ^ - 92 ) + (Ci92 - C 291 )) - 
8C?C9 - [(l-Ci)[9i-92^^-(9i-0(92-0] + C2[9i^^-(9i-0^]j‘ 

n _ +2 ^(12; 12)i-^^i 

9i [(Ci^ - 9i)^] [C 2 gi + Ci9i - (C 291 - Cl 92 )^] 

_^^_ 

1 - Cl [C 2 (Cl/ - 9 i)2 + Cl{C 2 l - 92)2 _ ^ _ ^^^^^^2 


X 


(37) 
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p2 ^ „4fCl^^Q2-21') C1C2 / [(1 - C2)qf + Clgl • g2][(l - Cl)g2 + C2gl ■ g2] 

’ (1 - Ci)(i - C2) [ [C2g? + Ci 92 - (C2<?i - Ci<?2)^]^ 

(Ci^ - qi) ■ [(Ci^ - 9 i) + (C291 - C192)] {C2I - 92) ■ [(C2^ - 92) + (Cl92 - C291)] 

^ [(Cl^ - 9 i) 2 ] [(C 2 ^ - 92 ) 2 ] 

91 • [(Ci^ - 91) + (C291 - Cl92)] {C2I - 92) • [(1 - Ci )92 + C291] 

9 l [(C2^ - 92)^] [C29? + C19I - (C291 - Ci 92 )^] [C2(Ci^ - 9 l)^ + Ci(C2^ - 92)^ - (C291 - Ci 92 )^] 

92 • [(C2^ — 92) + (Cl 92 — C291)] (Cl^ ~ 9 l) • [(1 — C 2 ) 9 l + Cl92] I 

9 i [(Ci^ - 9 i)^] [C29? + C19I - (C291 - Ci 92 )^] [C2(Ci^ - 9 l)^ + Ci(C2^ - 92)^ - (C291 - Ci 92 )^] J 

9 i • ((Ci^ - 9 i) + (C291 - C192)) 92 • (iC2l - 92) + (Ci 92 - C291)) 

9 i 92 ((C2^ — 92) + (Ci 92 — C291)) ((Ci^ ~ 9 i) + (C291 — C192)) 

(Ci^ - 9 i) • ((Ci^ - 9 i) + (C291 - Cl92)) {C2I - 92) • ({C2I - 92) + (Cl 92 - C291)) 

iCil - qi)'^{C,2l - 92)^((C2^ - 92) + (Ci 92 - C291)) ((Ci^ - 9 i) + (C291 - C192)) 

9 l • ((Cl^ - 9 i) + (C291 - C192)) {(,2l - 92) • ((C2^ - 92) + (Ci 92 - C291)) 

9 ?(C 2 ^ - 92)^ ((C2^ - 92) + (Ci 92 - C291)) ((Cl^ - 9 l) + (C291 - Cl92)) 

(Ci^ - 9 i) • ((Ci^ - 9 i) + (C291 - Cl92)) 92 ■ {{C2I - 92) + (Cl 92 - C291)) 

(Ci^ ~ 9 i )^92 ((C2^ — 92) + (Ci 92 — C291)) ((Ci^ ~ 9 i) + (C291 — Cl92)) 

r I 9i • ((Ci^ - 9i) + (C291 - C192)) 92 • (92 - (Ci 92 - C291)) 

9 ? 9 i [C29? + Ci 9 i - (C291 - Ci 92 ) 2 ] ((Ci^ - 9 i) + (C291 - C192)) 

9i • ((Ci^ - 9i) + (C291 - Cl92)) (C2^ - 92) • (92 - (Cl92 - C291)) 

9?(C2^ - 92)2 [C29? + C19I - (C291 - Ci 92 )^] ((Ci^ - 9 i) + (C291 - Cl92)) 

92 • ((C2^ - 92) + (Cl92 - C291)) (Ci^ - 91) • ((Ci^ - 91) + (C291 - Cl92)) 

92 (Ci^ - 9 i)^((Ci^ - 9 i) + (C291 - C192)) C2(Ci^ - 9 i)^ + Ci(C2^ - 92)^ - (C291 - Ci 92 )^ 

(C2^ - 92) ■ {{C2I - 92) + (Ci 92 - C291)) (Ci^ - 9 i) • ((Ci^ - 9 i) + (C291 - Cl92)) 

(C2^ - 92 )^(Ci^ - 9 l)^((Cl^ - 9 i) + (C291 - C192)) C2(Cl^ - 9 l)^ + Ci(C2^ - 92)^ - (C291 - Cl92)^ 


+ (Cl)9l ^ C2,92) 


(38) 



18 




C1C2 

f C1C2 

9 i 92 - (91 ■ 92)^ 


(1 -Cl)(l-C 2 )' 

, 9i9i [C29? + Ci9i - (C291 - Ci 92)^]^ 



1 

X 

to 

1 

to 

{(2I - 92) X (C192 - C291) 


[(Ci^ - 9 i)^] [(C 2 ^ - ^ 

2 ) 2 ] [C 2 (Ci?- 9 i )2 

9 i X (Ci^ — C192) 

+ Ci(C2^ - 92)2 - (C291 - Cl92)2 

{(21 - 92) X [(1 - Cl)92 + C291] 

2 


+ 


?! [(C2^ - 92 )^] [C29f + C 19 I - (C 291 - Ci92)^] [C2(Ci^ - 9l)^ + Cl(C2^ “ 92 )^ “ (C29l “ Cl92)^] 
92 X (^ 2 ^ — C 291 ) (Cl^ ~ 9l) X [(1 — C2)9l + Cl 92 ] 


9 i [(Ci^ - 9i)^] [C29? + C19I - (C291 - Ci 92 )^] [C2(Ci^ - 9l)^ + Cl(C2^ - 92)2 - (C291 - Ci 92 )^] 
9i X (Ci^ — C192) 92 X {C2I — C291) 


9?9 | ((C2^ - 92) + (Ci 92 - C291)) ((Ci^ - 9i) + (C291 - Cl92)) 

(Ci^ “ 9i) X (C291 — C192) {C2I — 92) X (Ci 92 — C291) 


(Cl^ - 9l)^(C2 ^ - 92)^ (^(C2^ - 92) + (Ci 92 - C291)) ((Cl^ - 9i) + (C291 - C192)) 
9i X (Ci^ — Cl92) {C2I — 92) X (Cl92 — C291) 


9 ?(C 2 Z- 92)2 

((C2^ - 92) + (Ci 92 - C291) 
(Ci^ - 91) X (C291 - Cl92) 

((Ci^- 9 i) + (< 
92 X (C2? - C29I 

(291 - Cl 92)) 

(Ci^ - 9 i 

?ql{{(,2i - 92) + « 
Cs 

192-C291)) ((Ci^ 
9 i X (Ci^ — Cl92) 

- 9 i) + (< 

92 X qi 

(291 - Cl92)) 

[i-Cil 

9 i 9 i [C29? + Ci 9 i - 

9 i X (Ci^ — Cl92) 

(C291 - Ci 92 ) 2 ] ((Ci^ - 9 i) + (C29 
{C2I — 92) X ^92 — (Ci 92 — C291)) 

- C192)) 


qUC2l - q2Y [C29? + C19I - (C291 - Ci 92)^] ((Ci^ - 91) + (C291 - Cl92)) 



92 X {C,2l - C291) 

(CiZ - 9 i) X (C291 - Cl92) 


9 i(Ci^ - 9 i) 2 ((Ci^ - 9 

l) + (C291 - Cl92)) 
(C2^ - 92) X (Ci 92 - 

C 2 (Cl? 

C291) 

- 9 i )2 +Cl(C2^ - 92)2 - ( 

(Ci^ - 9 i) ■ (C291 - C192) 

C291 - Ci 92)2 


{C2I - 92)^(Ci^ - 9i)^((Ci^ - 9i) + (C291 - C192)) C2(Ci^ - qiY + Ci(C2^ - 92)^ - (C291 - C192)" 
+ (Cl,9l ^ C2,92) 


Finally we have 


da 


d?kidk+d?k2dk+ 2 


-^(0 ^11 + ^11 + -f"l^2 + ^12 + (9lj Cl ^ 92, C2) 


(39) 


(40) 


Note that the partial symmetrization of some terms in Eqs. (I38I39I) should not be confused with the full symmetrization 
in the final expersion Eq. (l40l) . Indeed, there is no double counting. 
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A. The soft (eikonal) limit 


These expressions are fairly complicated, but it is straightforward to take the soft limit where the photons are much 
less energetic than the quark. This means C; Cii C 2 ^ 1- In this limit we have 


B{p-q,s,s',q,i) = g 



(41) 


Cip-qi - q 2 ,s,s',qi,i,q 2 ,j) = 


qiiq2j 


P+V C 1 C 2 + g?§) 


(42) 


These simple leading soft expressions give vanishing production amplitude. The reason is simply that they do not 
depend on the transverse quark momentum p, and therefore there is a complete cancellation in the amplitude Eq. O- 
Thus we proceed to expand to next order in We only need to keep the terms which depend on the momentum p, 
as other terms cancel just like with the leading eikonal pieces. Therefore, we obtain at next to leading order. 


B{p-q,s,s\q,i) = g^J^^ q^ - C 


^, p^-{p- q? 

q^-C\p^ - 5 -9* 


(5ss' = 


2 2 

Cp+ 


., 2 p • g 

qi-Q\Pt - —qi 


(43) 


C{p - qi - q 2 ,s,s',qi,i,q 2 ,j) = 


P+\ CiC2g2(g2 + g2|) 

f 2p-qi ^ ( 2 p-{qi+q 2 ) ^ 

quqij - Cl {p^ - —qii q 2 j - C2<?h Pj -5—57^^2^ 

V 9l / V 


Sss'. (44) 


In both of these expressions we have only kept terms proportional to the quark momentum p. To calculate the 
amplitude F we need 


[B{-qi,s,s”,qi,i) - B{1 - qi, s, s” ,qi,i)] B{1 - qi - 92 , s”, s', 92 , j) = 5^—r 


'Cl 1 


V C 2 qfqi 


k- 


21 ■ qi 


9? 


qii q 2 jS ss' • (45) 


C'(-<?i - q 2 ,s,s”;qi,i,q 2 ,j) - C{1 - qi - 92 , s, s”; gi, ( 72 , j) = 
1 


9 


P+V C 1 C 2 ql{ql + ql^) 
Thus in the current approximation we obtain 
F{ 0 ,l,qi,q 2 ,ij,s,s') = 


/- I 1 gi ^ 21 ■ {qi + q2) 

Cl ( -) 92j + C29ii Ij - 2 _i_ 2 C 2 92j 

9i / \ 92 + 'll ^ y 


(46) 


9^-t 


1 1 


P+ V C 1 C 2 9? 


Cl 


Cl f, qi ^ 

,?('* "‘V®’ + 


, 2 Z • gi j 

.2 _L * I * ^2 9ii ) 92i + 


91 


C 2 (, 2 Z • (gi + 52 ) 

9ii \ h 9 , 2 c, 92j 


„2 I „ 2 ^ 

92 + 9i 


„2 I „2ft 

92 + 9i 


^ss' • 

(47) 


To calculate the cross section we need to symmetrize this expression with respect to the external photon momenta and 
polarizations. The most interesting situation is when the two photons have very different rapidities. If the correlation 
function exhibits a genuine long range ridge-like structure, like in the di-hadron case, such a structure should survive 
in this limit. We therefore take zi <C Z 2 - In this limit we have 


F{0,l,ki,k2,i,j) + F{0,i,k2,ki,j,i) = 


9^-t 



■ ^1 

V -21Z2 



21-ki. 


Zl 


2Z2 


p2u2 


21 -ko 


^2 kli ) ^2j + + 1.21,2 \ 1.2 ^2^ ) ku + ^2 r,2 ( k 


Zl 


Ail ^2 


‘11 * {k\ ~\- /C2) 


kl 


ki, 


16 1 


^Zl<^Z2 


p+ Y Zl kfk2 


u- 


21 ■k2 




kii. 


(48) 
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It is now easy to calculate the production cross section by squaring Eq. (l48l) . 

da 5 ^ 16 ^ Z 2 1 

d?kidk+d'^k2dk+ ~ p+^ 

This obviously shows no angular correlations. Thus if the correlations exist, they can only be significant when and 
^2 are of similar magnitude. 

When calculating the cross section in p+p(A) collisions, the integration over Xq can pick a contribution from Z 2 
which is not very small. We should consider this possibility. Let us therefore relax our approximation and consider 
<C 1, but Z 2 not necessarily small. As a preamble to this we have 



C{p-qi - g 2 ,s,s',gi,i,g 2 ,i)ci<i = 



V '^2 9i g?(i-C2) 


qu [C2{p-qi)k - q2k] 


2-C2 
. C2 


dkjdss' 


• (50) 


C{p- qi- q 2 ,s,s',qi,i,q 2 ,j)c^<s:i = g — 


,44 Cl 


V C2 qj b 2 Ci(i - Cl) -{p- qi^Ci - 9 ?(i - Ci)] 


q2j [Cm - qu] 


2-Cl 


Cl 


dlidss' ss' 


(51) 


We now write the symmetrized expression that enters the amplitude 

C'(p qi q 2 j s, s , qi,i, q 2 , j)q^=ki , 92 =^ 2 ,Ci = 2:1 ,C 2 =Z 2 A C'(p q\ q 2 , s, s , qi, i, q 2 , j)q 2 =ki ,<?i=/c 2 ;C 2 =zi ,Ci =^2 


1 21<1 


44 /iT 1 

= -9 -T 


1 


p+ Y 2:2 kl kl{l - Z2) 

4 /i7 1 


kii[z2{p- ki)k - k2k\ 


2 — Z 2 


dkjdss' ^^kja g 

2 — Z2 


p+ Y Zl kl [p'^Z2{l - Z2) - (p- k 2 Yz 2 - fcf (1 - Z2)] 

1 


22 

ku [z 2 Pl - k 2 l] 


44 
= 9 -T 


p+ Y 2i kl [p2z2(l - 22) - (p- ^2)^22 - A:|(l - Z2)] 


ku [z 2 Pl - k 2 l] 


22 
2-22 


dljSss' S 


22 


dljSss' ^^Ij'^^ss 


(52) 


where the last equality follows since 21 <C 1 . 
Next consider 


D{p,l,s,s',qi,i,q 2 ,j) = B{p - qi,s,s'\qi,i) - B{p + l- qi,s,s'\qi,i) B{p + l- qi - 92 , s”, s', 92 , j)- 
We have 

D{p,l,s,s',qi,i,q 2 ,j)c^<g:i = 0{^/Cl) 


(53) 


D{p,l,s,s\qi,i,q2,j)c2^i = 9‘^^J^q2j 


2 -Cl 
Cl 

1 


dkidss' „/ 


Lp 2 Ci(l - Cl) - ip-qi)'^Ci - 9?(1 - Cl) 
1 


(p + 0 ^Ci(i - Cl) - (p + ^ - 9 i)^Ci - < 7?(1 - Cl) 


1 


(CiP/c - qik) 

(Ci(pfe + h) - qik) 


(54) 


Therefore 


D{p,l,s.,s ,qi,i,q 2 , j)qi=ki ,g 2 =fc 2 , 0 = 21 X 2=22 D^p,!, s, s , qi,i, q 2 j j)q 2 =ki ,gi=fe 2 ;C 2 = 2 i,Ci =22 

1 


A / 

^ Si 2, _3 

^kj Oss' '^^kj ^ ss' 

1 

T 2 

V 2l 

[ 22 J 

kf 1 


2 :i<l 

{z2Pk - k2k) 


1 


{p + 0 ^ 22(1 - 22 ) - (p + ^ - ^ 2)^22 - kl{l - 22 ) 


{Z2{pk + h) - k2k) 


(55) 


These expressions can be further simplified, since p = 0 in the collinear factorization approach. This is however not 
necessary. It is clear from the above expressions that in this limit the full amplitude F + only depends on the 
transverse momentum fci through the factor ki/kl- The cross section therefore again contains no correlations between 
kl and k 2 - This indicates that the correlations may appear only when both 21 and 22 are not too small. In Sec. Ill, 
we present our numerical results for the full expression Eq. (HOI) and show that they reproduce special cases considered 
in this Appendix. 
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B. Separating fragmentation and direct parts in a limiting case 


The cross-section of the prompt di-photon can be generally written as a sum of direct and fragmentation contri¬ 
butions. Note that for the general case, such an expression would be very lengthy (see Appendix A), and would not 
add any new insight, especially since here we stay away from the collinear kinematics. However, the general structure 
of collinear divergence in master equation ( 1221 ) . is similar to the simple case of the soft approximation considered 
in Ref. [s^. The cross-section for production of a quark with momentum q and two prompt photons with momenta 
ki and ^2 in the scattering of an on-shell quark with momentum p off a hadronic target, in a kinematic limit that 
1 ^ 1 , 2 ! < Ip — q\ (soft approximation), can be written in the following simple form in terms of fragmentation and direct 
parts [33|, 


^^Direct 


(Phd?'kiTdrii(P\i2Tdrj2 

_|_ kiT(l -|- 1/Ci) + ^2T\,Xg (jqr -I- kir/Cil) ) 

(^TT) Qt ^ 


[ 1 1 

(ktkt kt^ \ 

1 1 

^?T^ 2 T idrCl + CIIHT + kiT/Cl - k 2 T/C 2 P) 

{Okjj. Oklj.) 

l^tdrCi + + kir/Ci - k 2 T/C 2 p 


+ (^1 ^ Cl ^ C 2 )) (56) 


^^Fragmentation 




d'^hd'^kiTdriid'^k2Tdr]2 |k;iTC2 — k2TCiP 

X + 1 /Ci) + ^2T\,Xg (fciT/Ci) ) + {ki O k2,zi <->■ (2), 

^^Fragmentation 


1 P 

k^j, kl^O O 


Cl 


d?hd?'kiTdrj^-^^d?k2Td'q^^'^r^^W' 2(27r)^ 

X -I- 1/Ci) + k 2 r|, aig {kiT/zi) ^ -I- (fci o ^ 2 : ^1 ^ C 2 ); 

where Cq is the fractional electric charge of the projectile quark, and the function O is defined as 

O = kfrpk2 +k2Tk^, 


1 + (k+k+k^ ^ _ 

2 + 0{k+ei + ktCl) 


C 1 C 2 


-^71 72 /h{Cl,C2, dr) 


(57) 


(58) 


with ki = and k^ = In Eq. (1571) . and are the single photon and double photon 

fragmentation functions. The fragmentation parameters Ci and C 2 are defined in Eq. (1231) . The fragmentation part 
was separated from the direct di-photon cross-section by introducing a hard cutoff, the fragmentation scale pp, which 
defines the collinear singular part that is subsumed in the photon fragmentation contribution. In Eq. (EZl, the first 
and second term give the single and double photon fragmentation contribution, corresponding to the kinematics where 
only one photon or both photons is emitted almost collinearly with the outgoing quark. 
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